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ABSTRACT 

The inner region of the accretion disk around a magnetized neutron star is subjected to magnetic 
torques that induce warping and precession of the disk. These torques arise from interactions between 
the stellar field and the induced electric currents in the disk. We carry out a global analysis of warp- 
ing/precession modes in a viscous accretion disk, and show that under a wide range of conditions typical 
of accreting X-ray pulsars, the magnetic warping torque can overcome viscous damping and make the 
mode grow. The warping/precession modes are concentrated near the inner edge of the disk (at the 
magnetosphere-disk boundary), and can give rise to variabilities or quasi-periodic oscillations (QPOs) in 
the X-ray/UV/optical fluxes from X-ray pulsars. We examine the observed properties of mHz QPOs in 
several systems (such as 4U 1626-67), and suggest that some hitherto unexplained QPOs are the results 
of magnetically driven disk warping/precession. 

Subject headings: accretion, accretion disks - stars: neutron - stars: magnetic fields - pulsars - 
binaries: close 



1. INTRODUCTION 

Disk accretion onto a magnetic star occurs in a variety of 
astrophysical contexts, including accreting neutron stars, 
white dwarfs and pre-main-sequence stars (e.g., Frank et 
al. 1992). The basic picture of disk-magnetosphere in- 
teraction is well-known: at large radii the disk is unaf- 
fected by the stellar magnetic field; a somewhat sudden 
transition occurs when the stellar field disrupts the disk 
at the magnetospheric boundary (where the magnetic and 
plasma stresses balance), and channels the plasma onto the 
polar caps of the star. A large number of theoretical pa- 
pers have been written on the subject of the interaction be- 
tween accretion disks and magnetized stars (e.g., Pringle & 
Rees 1972; Ghosh & Lamb 1979, 1992; Aly 1980; Anzer & 
Borner 1980, 1983; Lipunov et al. 1981; Wang 1987, 1995; 
Aly & Kuijpers 1990; Spruit & Taam 1993; Shu et al. 1994; 
van Ballegooijen 1994; Lovelace et al. 1995,1998; Li, Wick- 
ramasinge & Riidiger 1996; Campbell 1997; Lai 1998; 
Terquem & Papaloizou 2000), but numerical study of this 
problem is still in its infancy (e.g., Hayashi et al. 1996; 
Miller & Stone 1997; Goodson et al. 1997; Fendt & Elstner 
2000). Outstanding issues remain, including the efficiency 
of field dissipation in/outside the disk, whether the disk 
excludes the stellar field by diamagnetic currents or the 
field can penetrate a large fraction of the disk, whether 
the threaded field remains closed (connecting the star and 
the disk) or becomes open by differential shearing, and 
whether/how magnetically driven wind is launched from 
the disk or the magnetosphere/corotation boundary. 

Many previous theoretical papers have, for simplicity, 
adopted the idealized geometry in which the magnetic 
axis, the spin axis and the disk angular momentum are 
aligned. However, in Lai (1999), it was shown that under 
quite general conditions, the stellar magnetic field can in- 
duce warping in the inner disk and make the disk precess 
around the spin axis (see §2). Such magnetically driven 
warping and precession open up new possibilities for the 



dynamical behaviors of disk accretion onto magnetic stars. 
In Shirakawa & Lai (2001) we have studied these effects 
in weakly magnetized accreting neutron stars and showed 
that the magnetic warping/precession effects may explain 
several observed features of low-frequency quasi-periodic 
oscillations in low-mass X-ray binaries. 

In this paper we study global, magnetically driven warp- 
ing/precession modes of inner disks of highly magnetized 
(B ~ 10 12 G) neutron stars (NSs), as in accreting X-ray 
pulsars (§2 and §3). Such study is the first step toward un- 
derstanding the observational manifestations of the mag- 
netic warping/precession effects. We are motivated by the 
observations of the milli-Hertz quasi-periodic oscillations 
(QPOs) in a number of X-ray pulsars (see Table 1 in §4 
and references therein). Of particular interest is the re- 
cent detection of the 1 mHz optical/UV oscillations in 4U 
1626-67 (Chakrabarty et al. 2001). In §4, we use the results 
of §3 and suggest that magnetically driven disk warping 
and precession naturally explain this and some other mHz 
variabilities observed in X-ray pulsars. 

2. MAGNETICALLY DRIVEN WARPING/PRECESSION AND 
ITS GLOBAL MODES 

The inner region of the accretion disk onto a rotating 
magnetized central star is subjected to magnetic torques 
that induce warping and precession of the disk (Lai 1999). 
These magnetic torques result from the interactions be- 
tween the accretion disk and the stellar magnetic field. 
Depending on how the disk responds to the stellar field, 
two different kinds of torque arise: (i) If the vertical stel- 
lar magnetic field B z penetrates the disk, it gets twisted 
by the disk rotation to produce an azimuthal field AB,p — 
T(B Z that has different signs above and below the disk 
(£ is the azimuthal pitch of the field line and depends on 
the dissipation in the disk), and a radial surface current 
K r results. The interaction between K r and the stellar 
B^ gives rise to a vertical force. While the mean force 
(averaging over the azimuthal direction) is zero, the un- 
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even distribution of the force induces a net warping torque 
which tends to misalign the angular momentum of the disk 
with the stellar spin axis, (ii) If the disk does not allow the 
vertical stellar field (e.g., the rapidly varying component 
of B z due to stellar rotation) to penetrate, an azimuthal 
screening current will be induced on the disk. This 
interacts with the radial magnetic field B r and produces 
a vertical force. The resulting precessional torque tends to 
drive the disk into retrograde precession around the stellar 
spin axis. 

In general, both the magnetic warping torque and the 
precessional torque are present. For small disk tilt angle 
(the angle between the disk normal and the spin axis), the 
precession angular frequency and warping rate at radius r 
are given by (see Lai 1999) 1 

"pW = -,.7^_yL^ *m (!) 



T w (r) = 



ir 2 r 7 fl(r)Y;(r)D(r) 

Cm 2 



cos 9, 



(2) 



47rr 7 f2(r)I](r) 

where \x is the stellar magnetic dipole moment, 9 is the 
angle between the magnetic dipole axis and the spin axis, 
O(r) is the orbital angular frequency, and £(r) is the sur- 
face density of the disk. The dimensionless function D(r) 
is given by 



D(r) = max ( ^/r 2 /rf n - 1, f2H(r)/r h 
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where H(r) is the half-thickness and r m is the inner radius 
of the disk. The function F(9) depends on the dielectric 
property of the disk. We can write 

F{6) = 2/ cos 2 (9 -sin 2 (9, (4) 

so that F(9) = — sin 2 9 if only the spin- variable vertical 
field is screened out by the disk (/ = 0), and F(9) = 
3 cos 2 9 — 1 if all vertical field is screened out (/ = 1). In 
reality, / lies between and 1. For concreteness, we shall 
set F{9) = — sin 2 9 in the following. 

For accretion-powered X-ray pulsars, the inner disk ra- 
dius ri n is given by the magnetosphere radius r m : 

1/7 



n 



GMM 2 



(3.4x 10 8 cm) 7] fi% 7 M~l /7 M~ 7 



-2/7 



(5) 

where M = (1AMq)Mia is the neutron star mass, 
M = (10 17 g s _1 )Mi7 is the mass accretion rate, /Z30 = 
/z/(10 30 Gem 3 ), and r\ ~ 0.5 — 1. For typical parameters, 
the precession frequency is 
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where we have used f2(r) = (GM/r 3 ) 1 / 2 , and r 8 = 
r/(10 8 cm). The warping rate T w (r) is of the same or- 
der of magnitude as Ct p (r). 

Since the precession rate fl p (r) depends strongly on r, 
coupling between different rings is needed to produce a 
global coherent precession. The coupling can be achieved 
either by viscous stress or through bending waves (e.g., 

1 Note that the stellar spin frequency f2 s does not appear in eqs. (1) 
averaged out; this is justified because Q s 2> \Q P \, \T W \. 



Papaloizou& Pringle 1983; Papaloizou & Tcrquem 1995; 
Larwood et al. 1996; Terquem 1998). In the viscosity 
dominated regime (i.e., the dimensionless viscosity param- 
eter a greater than H/r), the dynamics of the warps can 
be studied using the formalism of Papaloizou & Pringle 
(1983) (see also Pringle 1992; Ogilvie 1999; Ogilvie & 
Dubus 2001). Wc model the disk as a collection of rings 
which interact with each other via viscous stresses. Each 
ring at radius r has the unit normal vector l(r, t) . In 
the Cartesian coordinates, with the z-axis along the neu- 
tron star spin, we write 1 = (sin (3 cos 7, sin (3 sin 7, cos (5), 
with /3(r, t) the tilt angle and 7(7*, t) the twist angle. For 
f3 -C 1, the dynamical warp equation for 1 (Lai 1999; see 
Papaloizou & Pringle 1983; Pringle 1992) reduces to an 
equation for W(r, t) = (3(r, t)e l ^ r ^: 
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where J' = dj /dr, V\ is the usual viscosity , and v 2 is the 
viscosity which tends to reduce the disk tilt. We assume 
that the ratio of v 2 to v\ is constant. In deriving eq. (7), 
we have used the relations for the radial velocity and sur- 
face density: v r = -'iviJ~ x /2r and S = MJ/3irvi. The 
values and functional forms of v\, v 2 , Q p , T w [see eqs. (1) 
& (2)], and the dimensionless function J{r) [see eq. (16) 
below] depend on disk models. 

2.1. Power-law Disk Models 

To gain insight on the properties of the global warping- 
precessional modes, we consider power-law disk models, 
with £(r) oc r^. We also assume that D(r)= constant, 
3(f) = 1, ^2/^1 = 1, and Q(r) cx r _3//2 as in the Keple- 
rian flow. Then, from eqs. (1) & (2), wc have 
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Using the relation S = MJ ' /2>-kv\, the viscosity rate r vi g C 
can be written as 

-viscW = ^ = ^c(nn) ("f ) " 2 • (10) 

We look for a solution of the form W(r,t) — e %crt W(r) 
with the complex mode frequency o (= ay + ia^). It is 
convenient to define dimensionless quantities by 

0~ = fT(7 , in)Tvisc(^in), Q p = Q p (^*in)7visc (^in ) , 

= (r[ n )T v [ sc (ri n ) . (ff) 
Equation (7) then reduces to the dimentionless form: 
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where x = r/r- ln . 

& (2) since the variation of the field geometry due to the spin has been 
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It is clear from eq. (12) that, for a given /i, the mode fre- 
quency a depends only on two dimensionless parameters 
Cl p and f^. To solve eq. (12) for the complex eigenfunction 
W{x) and eigenvalue a, six real boundary conditions are 
needed. In our calculation, the disk extends from x; n = 1 
to x out = 50. For large x and large |<r|, equation (12) can 
be solved analytically, giving 
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where the sign of (ier) 1 / 2 should be chosen so that W(x) — > 
as x — > 00. This approximate analytical solution, evalu- 
ated at Xout, together with its derivative, gives four (real) 
outer boundary conditions. The inner boundary condition 
generally takes the form W'(xi n ) = aW(xi n ), with a be- 
ing a constant. Most of our results in §3 are based on 
a = (corresponding to zero torque at the inner edge of 
the disk), although we have experimented with different a's 
and found that for \a\ ^ 1 similar results are obtained (see 
Shirakawa & Lai 2001). In numerically searching a mode, 
we make a guess for the eigenvalue a and integrate eq. (12) 
from x out to x; n . Since W{x) changes very rapidly from 
x ou t to Xi n , we rewrite eq. (12) in terms of a new function w 
defined as W = e w and use that equation for integration. 
We find the correct value of a that satisfies the boundary 
conditions using the globally convergent Newton method 
(Press et al. 1992). 

2.2. Middle- Region Solution of the a-Disk 

Here we consider the "middle-region" (gas-pressure -and 
scattering-dominated) solution of the a-disk (Shakra & 
Sunyaev 1973; Novikov & Thorne 1973) which is relevant 
to the inner part of the disk in accretion-powered X-ray 
pulsars. In this model, 

£ = (7.5 x ltfgcm-^aZt'^Ml^Mlfrf^J^lU) 

H = (1.0 x lO'^aZl^M-J^M^rf^J 1 / 5 ,^) 

where a = (0.1)a_i is the a-viscosity parameter. The 
dimensionless function J{r) is given by 2 
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where £ is a dimensionless parameter with < £ < 1 
(£ = corresponds to to zero net angular momentum 
transfer across the inner disk, i.e., when the star is in spin 
equilibrium). Substituting equation (14) into equations 
(1) and (2), and using fi(r) = (GM/r 3 ) 1 / 2 , we get 
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the viscosity rate is calculated as 
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e lat W{r), and using the dimensionless 
and T w as defined in eq. (11), we can 
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write equation (7) in the dimensionless form: 
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where D in = D(r in ), J in = J(r in ), and J' = dj/dx. 
Using eqs. (3), (5), and (15), we can calculate An as 
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Using eqs. (5) & (21) in eqs. (17), (18), & (19), we can 
as 
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Note that under the assumptions of J7"(x) = 1 and D(x) = 
An, cq. (20) reduces to eq. (12) with ^ = —0.6. 

3. NUMERICAL RESULTS 
3.1. Mode Eigenfunction and Eigenvalue 

We first consider the power-law disk models of §2.1. For 
a given set of parameters (/x, Cl p , T w ), equation (12) allows 
for many eigenmodes. We shall focus on the "fundamen- 
tal" mode which is more concentrated near the inner edge 
of the disk and has larger oy (global precession frequency) 
and smaller Ui (damping rate) than any other "higher- 
order" modes. 

Figure 1 shows the tilt angle [3(x, t = 0) = |TU(x)| as- 
sociated with the modes for different sets of (Cl p , T w ), all 
with [i = —0.6 (corresponding to the "middle-region" a- 
disk with J{x) = 1 and D{x) = An)- We see that as 
\£l p \ (note Q p < due to retrograde precession) and T w 
increase, the modes become more concentrated near the 
inner radius of the disk. This behavior can be understood 
heuristically: for a given |f2 p (ri n )|, a larger \Q p \ implies 
smaller viscosity [see eq. (11)], and thus the coupling be- 
tween different disk radii is reduced, and the mode is less 
spread. 

Figure 2(a) shows the magnitude of the mode frequency 
(oy < due to retrograde precession) |5y| in units of \Q P \ 
as a function of |f2 p | for different values of T w . The ranges 
of |f2 p | and T w are chosen to be from 10 to 1000 to cover 
possible values of parameters for accretion-powered X-ray 
pulsars [see eqs. (22) & (23) with J- XVL — 1]. We include 
the (i = —1.0 and fj, = 1.0 results as well as the "middle 
region" /i = —0.6 result to show how our results vary with 
the change of the surface density power-law. We see that 
\a r /tl p \ = |ay/Qp(ri n )| always lies between 0.3 to 0.85. 
The ratio 
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2 Magnetic fields threading the disk can modify J{r) in a modcl-dcpcndcnt way (sec Lai 1999 for an example) 
still be approximated by cq. (16). 
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is consistent with the behavior of the mode eigenfunction 
(see Fig. 1) that a larger \Cl p \ and T w make the mode more 
concentrated near the inner disk edge. 

3.2. Global Warping Instability Criterion 

In the absence of magnetic warping torque (r^, = 0), we 
expect the disk warp to be damped by the viscous stress 
acting on the differential precession (f2 p ). This "mag- 
netic Bardeen-Petterson effect" (Lai 1999) is analogous 
to the usual Bardeen-Petterson effect (Bardeen & Petter- 
son 1975), where the combined effects of viscosity and dif- 
ferential Lensc-Thirring precession align the rotation axis 
of the inner disk with the spin axis of the rotating black 
hole (or rotating, non-magnetic NS). The competition be- 
tween the magnetically driven warping (T w ) and the mag- 
netic Bardeen-Petterson damping can be determined by 
our global analysis. 

Figure 2(b) shows the damping rate &i as a function of 
\£l p \ for different values of T w . We see that bi decreases 
as increases, and becomes negative (implying mode 
growth) when the ratio r„,/|f2 p | is sufficiently large. By 
solving eq. (12) with /x = —0.6, we find that the numerical 
value of di can be approximated by 

b % = -af w +b\n p \°- 6 , (24) 

with a ~ (0.5 — 1.0) and b ~ (0.5—1.0); this result is insen- 
sitive to modest change of the surface density power-law 
(/i = —1 to 1). For the mode to grow (a^ < 0) we require 

t w £ 2|ry ' 6 Global Warping Instability. (25) 

Using eqs. (22) and (23), this condition becomes 

xaZ^M^M^jr 034 ;> 1 (26) 

We see that for parameters which characterize X-ray pul- 
sars the mode growth condition can be satisfied, although 
not always. In general, high (but not unreasonable) C (> 
a few) and small J m (see §3.3) are preferred to obtain 
growing modes. 

3.3. Effect of J (r) 

The results of §3.1-3.2 are based on eq. (12), correspond- 
ing to eq. (20) with D(x) = D ln and J(x) = 1. Here we 
consider the solutions of eq. (20) with the function J(x) 
given by eq. (16). 

Figure 3 shows the mode frequency a as a function 
of Jin = 1 — £ [sec cq. (16)] for three different sets of 
(Q p , f w ): (-38JT 1 - 1 , 21JJ- 1 ), (-7.6JT 11 , 38JT 1 ), and 
(-38 jr 1 - 1 , 43J i ~ 1 ); these are obtained from cqs. (22) & 
(23) with (sin 2 (9, C)=(0-5, 5), (0.1, 5), and (0.5, 10), re- 
spectively, while other parameters being fixed to the stan- 
dard values [77 = 0.5, z^/^i = 1, M30 = 1, ot-i = 1, Mi. 4 = 
1, M17 = 1]. In the calculations, we set D(x) = An for 
simplicity; using the more accurate function D(x) given 
in eq. (3) only slightly changes the numerical results. We 
see from Fig. 3 that |<7 r /O p | is insensitive to the choice 
of £ since the most of the dependence on J7i„ is already 
absorbed into the definition of £l p [see eq. (22)]. We also 
see that a small J7i n tends to increase <7i/|0 p |, although 



growing warping modes still exist for a wide range of pa- 
rameters. We find that the simple global warping insta- 
bility criteria given in eqs. (25) and (26) can be used for 
Jin > 0.1. For smaller J m (<, 0.1), <7j should be obtained 
numerically to determine whether the mode grows or gets 
damped. 

4. APPLICATIONS TO MILLI-HERTZ QPO'S IN ACCRETING 
X-RAY PULSARS 

QPOs with frequencies 1 — 100 mHz have been detected 
in at least 11 accreting X-ray pulsars (Table 1; see also 
Boroson et al. 2000). These mHz QPOs are often inter- 
preted in terms of the beat frequency model (BFM; Al- 
par & Shaham 1985; Lamb et al. 1985) or the Keplerian 
frequency model (KFM; van der Klis et al. 1987). In the 
BFM, the observed QPO frequency represents the beat be- 
tween the Keplerian frequency vk at the inner disk radius 
n n and the NS spin frequency v s [i.e., ^qpo = ^(fm)— v s \- 
In the KFM, the QPOs arise from the modulation of the 
X-rays by some inhomogeneities in the inner disk at the 
Keplerian frequency [i.e., ^qpo = ^(Hn)]- However, we 
see from Table 1 that for several sources, more than one 
QPOs have been detected and the difference in the QPO 
frequencies is not equal to the spin frequency Thus it 
is evident that the KFM and/or the BFM cannot be the 
whole story. We also note that in both the KFM and the 
BFM, it is always postulated that the inner disk contain 
some blobs or inhomogeneities, whose physical origin is 
unclear. 

Here we suggest a "Magnetic Disk Precession Model" 
(MDPM) for the mHz variabilities and QPOs of accreting 
X-ray pulsars. The magnetically driven precession of the 
warped inner disk (outside but close to the magnetosphcrc 
boundary) can modulate X-ray/UV/optical flux in several 
ways: (i) The observed radiation (in UV and optical, de- 
pending on r; n ) from the inner disk due to intrinsic viscous 
dissipation varies as the angle between the local disk nor- 
mal vector and the linc-of-sight changes during the preces- 
sion; (ii) The flux of reprocessed UV/optical disk emission 
visible along our sight also varies as the reprocessing ge- 
ometry evolves; (iii) Modulation of the X-ray flux arises 
from regular occulation/obscuration of the radiation from 
the central NS and magnetosphcrc by the processing in- 
ner disk. In the MDPM, we identify fQpo with the global 
precession frequency driven by the magnetic torques. Our 
calculations in §3 show that under a wide range of condi- 
tions, the warping/precession mode is concentrated near 
the disk inner edge, and the global mode frequency is equal 
to A = 0.3 — 0.85 (depending on details of the disk struc- 
ture) times the magnetically driven precession frequency 
at nn = r m . Thus we write ^qpo = -A|fip(>"in)|/27r. Using 
cq. (17) together with eqs. (5), (14), and (21), we have (for 
the a-disk model), 

^qpo = (15.7 mHz)A sin 2 9 i4 al*M^l /W M^ /5 
( r in r 49/io/ An N-i 
ll0 8 cnJ \0.1 J Jin 

/ 77 \ -4.93 

= (0.83 mHz) A (^J sin 2 9 

xa™ 5 M°-j 8 Mb n jr - 7 . (27) 
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We also note that the Keplerian frequency at r — r nl = r rn 
is 

Mnn) - (985 mHz) (^) ~%^ 6/ VM/ 7 . (28) 

We see from eq. (27) that the MDPM can produce QPOs 
with frequencies of order 1 mHz; larger values of z^qpo 
would require J ln <C 1 (corresponding to low surface den- 
sity at r; n ). The value of J- m depends on details of the 
physics at the inner edge of the disk, therefore is uncer- 
tain. Let Vr(r; n ) = x c s{i"m), where V r — M/(2nrT,) is the 
radial velocity and c s = HQ is the sound speed of the disk, 
we find 




(29) 

Setting x = 1 would give the minimum value of Ji n . 
We now discuss several individual sources. 



4.1. 4U 1626+61 

The LMXB 4U 1626+67 consists of a v s = 130 mHz X- 
ray pulsar in a 42 min orbit with very low-mass (<J 0.1 M Q ) 
companion (see Chakrabarty 1998). QPOs at 48 mHz 
(and oscillations at 130 mHz) have been detected simul- 
taneously in X-rays and in the optical/UV band (Shinoda 
et al. 1990; Chakrabarty 1998; Chakrabarty et al. 2001). 
Thus it is natural to attribute the optical/UV variability 
to the reprocessing of the variable X-ray emission in the 
accretion disk. Since v s > 48 mHz, the KFM is problem- 
atic because the propeller effect would inhibit accretion 
when v s > (r m ). (It is still possible to ascribe the QPO 
to Keplerian motion at some radius farther out in the disk, 
but this is rather ad hoc.) The BFM is a viable alternative 
for the 48 mHz QPO. 

Recent HST observations by Chakrabarty et al. (2001) 
revealed a strong QPO around 1 mHz (centroid frequency 
in the range of 0.3-1.2 mHz, and Q = v/Ai/ of order 10) in 
the optical/UV band. This QPO is absent in simultaneous 
X-ray data, and is stronger in UV and weaker in the opti- 
cal band. These features can be naturally explained as due 
to warping of the inner accretion disk (see Chakrabarty 
et al. 2001). Indeed, using B ~ 3 x 10 12 G (from Bep- 
poSAX observations of a cyclotron feature; Orlandini et 
al. 1998) and Myj ~ 1 (corresponding to X-ray luminos- 
ity L x ^ 10 37 erg/s; Chakrabarty 1998), eq. (27) yields 
^qpo 0.34A(r//0.5)- 4 - 9 sin 2 6»a°;f JT ' 7 mHz, which 
can easily give the observed j^qpo = 1 mHz provided we 
allow for J7in < 1 [sec eq. (29)]. If we interpret the 48 mHz 
QPO with the BFM, we have vk{tui) = 178 mHz and 
r\ ~ 0.83, we therefore require J m <, 0.01 (depending on a 
and sin 2 9) . 

4.2. Other Sources 

For the other sources listed in Table 1, no "smoking- 
gun" signature of warped disk has been observed. How- 
ever, the numerical values of QPO frequencies indicate 
that the MDPM may be at work. 

Her X~l: This well-studied binary X-ray pulsar {y s = 
808 mHz) shows QPOs in the UV and optical bands at fre- 
quencies of 8 ± 2 and 43 ± 2 mHz; these QPOs most likely 
arise from the reprocessing of the disk oscillations by the 
companion star (Boroson et al. 2000). A QPO at 12 mHz 



in X-rays is also present, but its connection with the 8 mHz 
QPO is not clear (Moon & Eikenberry 2001b). Since 
v s > vqpo, the KFM is not applicable. The BFM predicts 
vk (fin) = 816 mHZ and 851 mHZ for z^qpo = 8 mHz and 
43 mHz respectively. For ^30 ~ 3 (see Table 1) , Mi 7 ~ 1 
(corresponding to Lx — 10 37 ergs; Choi et al. 1994), and 
M1.4 ~ 1, we have v K (r- m ) ~ 400 (77/O.5)- 3 / 2 mHz. Thus 
to explain the 8 mHz or 43 mHz QPO with the BFM 
requires 77 < 0.5. The phenomenology of the 12 mHz X- 
ray QPO is consistent with disk precession (see Moon & 
Eikenberry 2001b); this may be explained by the MDPM. 

LMC X~4: This persistent X-ray pulsar (y s — 74 mHz) 
exhibits large X-ray flares. QPOs at frequencies of 0.65- 
1.35 mHz and 2-20 mHz have been found during such flares 
(Moon & Eikenberry 2001a). Since v s > z/qpo, the KFM is 
not applicable. The BFM predicts ^(r m ) ~ 75 mHz and 
~ 76-94 mHz for z^qpo = 0.65- 1.35 mHz and 2-20 mHz, 
respectively. Equation (28) with the measured ^30 — 11 
(see Table 1) and M17 ~ 25 (corresponding to Lx — 
5 x 10 38 ergs; Moon & Eikenberry 2001c) gives !/^(r in ) = 
502 (77/0.5)^ 3 / 2 mHz, which is much larger than the values 
of vk{j 'in) required by the BFM (even for r\ = 1). Hence, 
it is difficult, if not impossible, to identify the observed 
QPO frequencies (especially for z^qpo = 0.65— 1.35 mHz) 
with the beat frequency. On the other hand, eq. (27) gives 
z/ QPO 1.2 A(r]/0.5)- 4 - 9 sin 2 ea°:l 5 jr - 7 mHz, which is 
close to the observed vqpo = 0.65 — 1.35 mHz. 

4U 0115+63: This transient source shows a broad 
62 mHz QPO during a flaring state (Soong & Swank 1989). 
Recent XTE observation also reveals a prominent QPO in 
the X-ray flux at 2 mHz (Heindl et al. 1999). Heindl et 
al. (1999) noted that this 2 mHz QPO may be explained 
by occultation of the radiation beam by the accretion disk. 
Analogous to the 1 mHz QPO of 4U 1626+67 (see §4.1), we 
suggest that magnetically driven disk warping/precession 
may be responsible for such occultation. With the well- 
measured magnetic field (B ~ 1.3 x 10 12 G) of the NS 
from the cyclotron lines, similar constraints on the system 
parameters can be obtained. 

For other sources listed in Table 1, either only sin- 
gle QPO is known, or there is no independent measure- 
ment/constraint on /U30 and M, thus the theoretical inter- 
pretation is currently ambiguous. 

5. CONCLUSIONS 

We have shown in this paper that the inner region of 
the disk around a strongly magnetized (~ 10 12 G) neu- 
tron star can be warped and will precess around the stellar 
spin (see Lai 1999). These effects arise from the interac- 
tions between the stellar field and the induced currents 
in the disk (before it is disrupted at the magnetosphcrc 
boundary). We have carried out a global analysis of the 
warping/precession modes and found that growing modes 
exist for a wide range of parameters typical of accreting 
X-ray pulsars. We therefore expect that the magnetically 
driven warping/precession effect will give rise to variabili- 
ties/QPOs in the X-ray/UV/optical fluxes. We have sug- 
gested that some mHz QPOs observed in several systems 
(e.g., 4U 1626-67) are the results of these new magnetic 
effects. Although there are significant uncertainties in the 
physical conditions of the magnetosphere-disk boundary 
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[and these uncertainties prohibit accurate calculation of 
the QPO frequency; see eq. (27)], we emphasize that the 
existence of these effects is robust. Continued study of the 
variabilities of X-ray pulsar systems would provide useful 
constraints on the magnetosphere-disk interactions. 
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Table 1 

Accretion-powered X-ray pulsars with mHz QPOs 



System 


spin frequency 


QPO frequency 


magnetic moment 


References a 




v a [mlizj 


^QPO [mHzJ 


,, b 
^30 




4U 1626-67 


130 


1, 48 


3.3 


[1], [2], [3], [4]t 


Her X-l 


807.9 


8, 12, 43 


3 


[5], [6], [7]t 




i zLi n 


(OUJ 






Cen X-3 


207 


35 




[9] 


LMC X-4 


74 


0.65-1.35, 2-20 


11 


[10], [11]+ 


4U 1907+09 


2.27 


55 


2.5 


[12], [13]+ 


EXO 2030+375 


24 


187-213 




[14] 


4U 0115+63 


277 


2, 62 


1.5 


[15], [16] 


XTE J1858+034 


4.5 


111 




[17] 


V 0332+53 


229 


51 


2.5 


[18], [19]+ 


A 0535+262 


9.71 


27-72 


9.5 


[20], [21]+ 



"References: [1] Shinoda et al. 1990; [2] Chakrabarty 1998; [3] Chakrabarty et al. 2001; [4]t Orlandini et al. 1998; [5] Boroson et al. 2000; [6] 
Moon & Eikenberry 2001; [7] + ; Makishima et al. 1999; [8] Wojdowski et al. 1998; [9] Takeshima et al. 1991; [10] Moon & Eikenberry 2001; 
[11]+ La Barbera et al. 2001; [12] in't Zand et al. 1998; [13] 1 " Makishima & Mihara 1992; [14] Angelini et al. 1989; [15] Soong & Swank 1989; 
[16] Hcindl et al. 1999; [17] Paul & Rao 1998; [18] Takeshima et al. 1994; [19] t Makishima et al. 1990; [20] Finger et al. 1996; [21] t Grove et 
al. 1995. These references include QPO discovery (no mark) and magnetic field strength B estimated from cyclotron features (with f). 

'The magnetic moment is calculated by /j, = BR 3 , or ^30 = Bi2R§. We assume R = 10 km, or Rq = 1 here. 
c marginal detection 
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Fig. 1. — The disk tilt angle of the warping/precession modes as determined from eq. (12) with fi = —0.6. The curves represent 
the fundamental modes for (Op, T m ) =(—10,100), (—100,10), and (—10,10) from left to right, with the corresponding mode frequency 
<7 = (oy,cr;) = (—6.9, —57), (—62, 12), and (—4.2, —1.7). The eigenfunction is normalized such that the maximum tilt angle is 1. 
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Fig. 2.— The upper panel shows the magnitude of the mode frequency |t7y | in units of as a function of for different values of t*w 
Different surface density power-laws, £ oc r M with /i = — 1, —0.6 and 1, are adopted [see eq. (12)]. The lower panel shows the corresponding 
mode damping rate &i (in units of |f2 p |). Note that negative cr; implies growing mode. 
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0.8 
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Oh 



<b 1.5 



-2 
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0.8 



Jin=W 



Fig. 3. — The upper panel shows the magnitude of the mode frequency oy| in units of |f! p | as a function of J\ 
sets A, B, and C: (Cl p , f w )= (-38JV 
panel shows the corresponding mode damping rate &i (in units of \O p \). Note that negative <jj implies growing mode. 



1 — £ for diffcrct parameter 



11 , HJin 1 ), (-t-GJfo 1 - 1 , 38 jr 1 ), and (-38 jr 11 , 43 1 7 in 1 ), respectively (see the text). The lower 



